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SUMMARY 


Reverse -flow theorems in aerodynamics are shown to be based on the 
same general concepts involved in many reciprocity theorems in the physi 
cal sciences. Reciprocal theorems for both steady and unsteady motion 
are found as a logical consequence of this approach. No restrictions on 
wing plan form or flight Mach number are made beyond those required in 
linearized compressible-flow analysis. A number of examples are listed, 
including general integral theorems for lifting, rolling, and pitching 
wings and for wings in nonuniform downwash fields. Correspondence is 
also established between the build-up of circulation with time of a wing 
starting impulsively from rest and the build-up of lift of the same wing 
moving in the reverse direction into a sharp-edged gust. 

INTRODUCTION 


Some of the most important results in the recent study of wing 
theory have been achieved through the development of reverse-flow rela- 
tions. The theorems already obtained are of outstanding practical util- 
ity and it appears obvious that the fullest exploitation of the methods 
has yet to be accomplished, either from a purely theoretical standpoint 
or in the routine calculation of wing characteristics. Attention to such 
problems in aerodynamics was initiated by von Karman (reference 1) who 
first announced the invariance of drag with forward and reversed direc- 
tions of flight for a nonlifting symmetrical wing at supersonic speed. 
Subsequently, advances in the theory were made by Munk, Hayes, Brown, 
Harmon, and Flax (references 2 through 7). Up to the present time, the 
most general results have been expressed by Ursell and Ward (reference 8') 
and by R. T. Jones (reference 9) in his attack on the study of wing 
shapes of minimum drag. In the forms given in the two latter papers the 
derived equivalences could be termed reciprocal or reciprocity relations 
rather than reverse-flow relations; in fact, this change in terminology 
divorces attention, momentarily, from the purely aerodynamic aspects of 
the results and, in this way, suggests a reorientation in terms of the 
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various similar relations appearing in other engineering fields. In the 
theory of elasticity, for example, a reciprocity theorem for small dis- 
placements of an elastic medium is so expressed as to appear in formal 
agreement with the statement of the result, given by Ursell and Ward (see, 
e.g., reference 10). This theorem is attributed to E. Betti and was pub- 
lished in 1872. A generalization was given by Lord- Rayleigh in 1873, and 
in various sections of his two volumes on the theory of sound (refer- 
ence 11) discussions of reciprocal relations in an elastic medium and for 
acoustic sources are given. In 1886, von Helmholtz (reference 12) 
obtained by means of variational methods applied to Hamilton’s character- 
istic function a reciprocal theorem for small changes in the momenta and 
coordinates of a general dynamical system in forward and- reversed motion. 
This result was commented on, in turn, by Lamb (reference 13) and an 
independent proof based upon Lagrange's generalized equations was given. 
The paper by Lamb is of particular interest since it contains the essen- 
tial idea underlying the development of reverse-flow theorems in wing 
theory. Thus, Lamb remarks, as had Lord Rayleigh previously, that reci- 
procity relations between sound sources do not apply directly in a moving 
atmosphere. He points out, however, that the reciprocity can be restored 
if the direction of the wind is also reversed. Further examples of recip- 
rocal theorems appear in electrostatics and in optics. 

The generality in the statement of reciprocity relations appears, . 
almost universally, to have held back their application to problems for 
which they are obviously, in retrospect, particularly fitted. This 
generality is even more apparent in some of the conclusions of Lord 
Rayleigh and von Helmholtz which apply to nonconservative systems . 

The purpose of the present paper is twofold. First, a close con- 
nection will be established between reverse-flow theorems in subsonic and 
supersonic, steady-state wing theory and known reciprocity relations 
between two solutions of the equation governing the flow field. In this 
way, machinery will be provided whereby extensions of existing results to 
the case of unsteady motion follow directly. Second, a number of partic- 
ular problems in wing theory in steady and unsteady flow will be consid- 
ered. It will be shown that, provided attention is limited to force and 
moment characteristics, the complexity of many solutions involving non- 
uniform flow fields, control surface deflections, and unsteady motion can 
be reduced considerably. In some cases, previously obtained solutions 
will be calculated. Comparison with the original calculations will 
almost invariably highlight the economy of effort in obtaining the final 
result. The utility of reverse-flow theorems is based on the fact that 
they build from known solutions and thus avoid the necessity of starting 
each problem anew. 
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GENERAL ANALYSIS 

RECIPROCITY RELATIONS FOR A CLASS OF PARTIAL 
DIFFERENTIAL EQUATIONS 


In this section, integral relations associated with linear partial 
differential equations will he reviewed from the standpoint of relating 
independent solutions. The subject matter is precisely Green's theorem 
and, in common with the usual expression of the theorem, it is preferable 
to treat the variables initially as abstract quantities. Consider, 
therefore, a class of linear partial differential equations of Becond 
order with independent variables Xj, X 2 , X m that may be thought of 

as rectangular coordinates in a space of m dimensions. Denote differ- 
entiation of the function i|r(Xi, X 2 , •••> %) with respect to the variables 
X^ and Xj by the subscript notation 


(♦) 


_ 

i j. dXjdXj 


and consider the differential equation 

m m 

L (t) = AijCtJij + B* = 0 


( 1 ) 


where, for the purposes at hand, Ajj = Aji are made independent of 

Xj and Xj, and B is a constant. Such equations fall within the class 
of self-adjoint equations. 


By Green's theorem (see, e.g., reference 14), it is possible to 
relate two arbitrary functions t and G by the integral expression 


P J [flL(t) - ♦!(«)] dV = - J J[ fiD n t - tonQldS (2) 


where the left member is a volume integral over a prescribed region in 
m-dimensional space and the right member is a surface integral over the 
hypersurface S enclosing the given region. Equation (2) certainly 
holds for any region in which ^ and Q and their first and second 
derivatives are continuous. The directional derivatives D n are defined 
in terms of the direction cosines ni, n^, ...,:% of the normal to the 
surface S with the stipulation that the normal is directed into the 
given region. Thus, setting 
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where v x , Vz, . . . , v m are the direction cosines of a line termed the 
"conormal," the directional derivative is defined by the expression 

Bn* = N $ V^lOi = N(i) v = N — (3) 

i=l £v 


If, finally, ^ and ft are assumed to satisfy equation (l), the left 
side of equation (2) vanishes and the resulting expression 

BD n \|rdS = J J ^D n fldS (4) 

is a general reciprocity relation expressing the functional dependence 
between two arbitrary solutions of equation (l). 

An interesting interpretation of equation (4) has been given by 
several writers (see, e.g., reference 15) and applies to the particular 
case when t and ft are identified with the perturbation velocity poten- 
tial Cp in the theory of incompressible-fluid flow, 
equation is Laplace's equation in three dimensions 

%cx + ^yy + ^zz = ® 

and the reciprocity relation takes the form, for two 
9 and 9’, 

ffozt dS °// <p 's dS 

where the directional derivatives are now along the true normals to the 
surface S enclosing a three-dimensional volume. It is known that any 
actual state of motion of a liquid for which a single-valued velocity 
potential exists can be produced instantaneously by the application of a 
properly chosen system of impulsive pressures. These impulsive pressures 
are directly proportional to the velocity potential^ plus an arbitrary 
constant which may in the present case be associated with the pressure 
of the uniform stream. Equation (6) is thus seen to represent summations 
of the cross products of impulse and normal velocity in two possible 
motions of a conservative system and is a special case of the dynamical 
theorem (reference 11) 


The governing 

(5) 

possible solutions 

( 6 ) 



00 • f. 00 1 • / r—j \ 

Z P r <lr = Z Pr <lr (7) 

r=i r=i 

where p r , 4r an d Pr* > 4r' are generalized components of impulse and 
velocity in any two possible motions of a system. 

The interpretation of equation (6) that leads to equation (7) pro- 
vides an indication of the close connection between reciprocal theorems 
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based upon the principles of least action and the symmetric character of 
Green's theorem for certain second-order differential equations. In the 
subsequent applications it will be convenient to proceed directly from 
equation ( 4 ) and seek to establish reciprocal relations between flow 
fields in wing theory. Such a process is well known when t or H is 
replaced by the elementary solution associated with a unit source and, 
in this case, establishes a general solution in terms of source and 
doublet distributions determined by arbitrary boundary conditions. The 
present objective is, however, different in that one wishes to get a 
symmetrical dependence between two general solutions. Moreover, the 
apparent symmetry of equation ( 4 ) must be consistent with physical con- 
siderations so that, for example, the velocity at point A induced by a 
single source at point B is equal to the velocity at point B induced 
by a source at point A. If, in the two systems, the effective free 
streams or flight directions are opposed, a fore and aft symmetry occurs 
•and the possibility of maintaining symmetry in the reciprocal relations 
becomes feasible. . . 


RECIPROCAL THEOREMS IN WING THEORY 
Reverse Flow for Subsonic Wings 

Consider a thin wing at possibly a small angle of attack and situ- 
ated in the immediate vicinity of a flat surface which is designated the 
xy plane. For sufficiently small thickness and angle of attack the 
perturbation velocity potential or any of the perturbation velocity com- 
ponents of the wing satisfy the linearized partial differential equation 
of compressible flow „ 

P %cx + ^yy + ^zz = ® ' (®) 

where |3 2 = 1 -Mq 2 = l-(U 0 /a 0 ) . Equation (8) applies in forward or reverse- 
flow provided the corresponding free -stream Mach numbers Mr, and M n ’ are 
equal. . In the accompanying sketch a 
lifting wing with plan form P is 
indicated along with the vortex wakes 
as the configuration would appear if 
the two flow fields were superimposed. 

It will be assumed that the wing chord 
is finite and that the profiles are 
closed. 

It is proposed to apply the gen- 
eral reciprocity relation to these 
flow fields in a manner similar to 
the approach used in developing the 
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basic solutions of the differential equation (see, e.g., references 16 
and 17). Thus, for subsonic flow, hemispherical regions of large radius 
and lying first above and then below the plane of the wing will be chosen 
as the volumes of integration. The surface integrals will therefore 
extend over a hemisphere with center at the wing and a flat surface that 
lies immediately adjacent to the z = 0 plane. This latter surface is 
subsequently to be brought into coincidence with the plane of the wing 
but must be considered first in its displaced position since only then 
can the flow be assumed free of possible singularities in perturbation 
velocities and their gradients. As in wing theory, in general, the 
attenuation of the perturbation potential and its gradients may be 
assumed of such a nature that the integrated contributions • of the wing 
and its wake over the hemispherical surfaces vanish in the limit as the 
radius becomes infinitely large. It remains therefore to consider the 
integrals over the surfaces at z = 0 + € and z = 0-e. Denote these 
surfaces as cr u and aj, respectively, where the subscripts u and l 
specify values above and below the z = 0 plane. Equation (4) can then 
be written as 



n directed upward 



where the primes denote conditions for 


n directed downward 
the reversed flow. 


(9a) 

(9b) 


In equations (9), ^ is now replaced by the perturbation velocity 
potential 9(x,y, z) and by the x or stream-wise component of per- 
turbation velocity u(x,y,z). * By virtue of the irrotationality of the 
flow, the gradients of u and w are related by the expression 

Bu _ + <5w 
Bn Bx 


and equations (9) can both be written in the form 



Integration by parts over either of the surfaces a u or a? gives 




( 10 ) 


At x = -oo the potential cp for the forward flow vanishes and at 
x = +“ the upwash w' in the reversed flow vanishes so that the first 
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term on the left is zero. The remaining double integrals have for their 
surfaces of integration the displaced planes o u and <J Z . In order to 
obtain a concise form of the reverse -flow theorem, it suffices to sub- 
tract the integrals extended over a z from the integrals over a u and 
let the planes approach coincidence with the plane of the wing. Since 
w,w' and u,u* are continuous everywhere except possibly in the immediate 
vicinity of the wing, the integration region can be restricted to planes 
slightly above and below the wing but extending beyond the wing edges. 
Provided the singularities at the edges can be disregarded, the analyti- 
cal expression of the steady-state reverse -flow theorem of Ursell and 
Ward (reference 8) becomes 



u z w z ')dS 



u 2 *w z )dS 


(ID 


for either lifting surfaces or symmetrical wings where P is the plan 
form of the wing in the z = 0 plane. 


It remains to discuss the effect of edge singularities. In the case 
of a lifting surface., square-root singularities in both u and w can 
occur at the leading and trailing edges -just on and off the wing, 
respectively. In the combined flow fields, the limiting process would 
then yield residue terms analogous to the leading-edge thrust of a lift- 
ing plate. If, however, the Kutta condition is imposed at the trailing 
edge for both the flow in forward and reverse direction,, a combination of 
singularities does not occur and equation (ll) is valid. If the leading 
edge of a symmetrical wing is rounded so as to produce a square -root 
singularity in w on the wing, a square-root singularity in u occurs 
just off the wing and a term corresponding to the leading-edge drag 
(reference l8) appears. If the geometry of the wing is fixed in the for- 
ward and reverse flow, however, the effect of these terms is canceled and 
equation (ll) still applies. It is important to point out that if the 
Kutta condition is not satisfied then the area of integration of 
equation (ll) cannot be restricted to the plan form. 


The two sides of equation (ll) are expressed in terms of the same 
coordinate system but it is usually preferable to associate with each of 
the two streams an x axis extending in the stream direction. To this 
end introduce now the subscripts 1 and 2 to denote the forward and 
reverse flow and the two coordinate systems. Thus, in general. 


yi = -y 2 + z i = z 2 


X 1 = “ X 2 + £> 


( 12 ) 
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where % and r| are arbitrary constants, and equation (ll) then becomes 

ff «u 1 (xi»y 1 )wu^* 1 >y 1 ) - u z x ( x i>yi) w z 2 ( x i^yi) J dy x = 

Pi 

u u 2 ( x a>y2) w Ul ( x 2>y2) - u 2 2 ( x 2>y2) w z 1 ( x 2^y 2 ) J dx 2 d y 2 (13) 



In the case of a symmetrical nonlifting wing, theory yields the 
relations 


"u = U Z> w u = “ w 2 


and for lifting surfaces 

u u = -Uj, 


w u = w 2 


It follows that in either case, equation (13) reduces to the form 




u 2 ( x 2 ,y 2 )Wi( x 2 ,y 2 ) dx 2 dy 2 (lU) 


where the velocity components can be evaluated on either the upper or 
lower surface of the plan form. Furthermore, the linearized pressure 
relation is 


P-Po = -Po u o u (15) 

where p is local static pressure and p Q is static pressure in the 
free stream. Since the wing profiles are assumed closed, equation (l^) 
becomes 

>yi) w 2 ( x i>yi) <ixi<ayi = 

p 2 

If, instead of specifying boundary conditions in a single plane, it 
is necessary to treat boundary conditions for a system of planes, the 
expression of the reverse-flow theorem is of the same general form as 
equation (ll) . Provided the Kutta condition is imposed at the trailing 
edges of all lifting surfaces, the relation becomes 


ff* 
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where the area of integration E extends over both, sides of all the 
wing surfaces, V n is the component of perturbation velocity normal to 
and directed away from each wing, and the subscripts 1 and 2 refer to 
forward and reverse flow in the two axial systems. 


Reverse Flow for Supersonic Wings 


The development of a reverse-flow theorem for supersonic wings 
parallels closely the analysis for the subsonic case. For either planar 
or multiplanar systems the conormal in equation (^) is, in fact, the nor- 
mal so long as the surface of integration is a plane parallel to the 
x axis. In the case of the single wing, for example, equations (9) 
apply where the surfaces a are slightly removed from the plane of the 
wing and where t satisfies the differential equation 

(Mo 2 -l) - tyy - t ZZ = 0 

In the limit as o approaches the z = 0 plane, the reversibility 
theorem takes the form of equations (ll) and (l6) provided the integra- 
tion extends beyond the edges of the wing- It is necessary to include 
these edges for wings with subsonic leading and trailing edges since 
singularities occur in the perturbation components and the solutions are 
not necessarily unique. For supersonic type edges, the area of integra- 
tion can be confined to the plan form of the wing and this is also true 
for subsonic edges provided the Kutta condition holds for all subsonic 
trailing edges in both the forward and reverse flow. Equation (l7) 
relates the two possible flows in the case of multiplanar systems. 


Reverse Unsteady Motion 


In the case of unsteady motion at either subsonic or supersonic 

flight speed, the basic equation may be taken in the form 

✓ 

ttt - *xx - +yy - +zz = 0 ( l8 ) 

where t = a 0 t ! , a Q is the speed of sound in the undisturbed region of 
the field, t\ is time, and \|/ is the perturbation velocity potential 
or any of the perturbation velocity components. Equation (l8) is the 
acoustic equation for small disturbances in three space dimensions and 
holds for a system of Cartesian coordinates fixed relative to the undis- 
turbed air. In applications to wing theory, therefore, * the wings move 
relative to fixed axes. 



10 


MCA TN 2700 


In the derivation of a useful theorem it is convenient to treat thin 
wings at small angles of attack and to assume that the motion takes place 
in the xy plane. The visualization of the time and geometry relations 
is relatively easy for two-dimensional wings moving at a uniform speed, 
as indicated in the sketch. The airfoil starts at time t' =0 and 

moves to the left at a constant 
velocity U 0 so that the trace of 
* the leading edge in the xt plane 
is x = -U 0 t' = -Mot and the 
trailing -edge trace is x = 2a - Mot. 
The lines x = ±t and x = 2a ±t are 
the traces of the extremities of the 
regions affected by the acoustic 
waves set in motion at t' = 0 by 
the leading and trailing edges. 

In the sketch, the wing has traveled 
a time T' = T/a 0 and the boundary 
condition determining the wing shape 
during the motion will be fixed by 
prescribing the value of vertical induced velocity w over the region 
"swept-out" of the xt plane by the wing. In order to determine a 
reciprocal theorem, a second wing is assumed to start at the final posi- 
tion of the first wing and to move with negative velocity until it has 
reached the initial position of the first wing. With these concepts in 
mind, it follows from equation (4) that the relation 

//"’ IS to4t ‘ If If wdxdt (19) 

holds where 9 and cp * are, respectively, the perturbation potentials 
for the forward and reverse motions. The region of integration is deter- 
mined by the area occupied in the xt plane by the wing, and the Kutta 
condition is assumed to apply to the trailing edge when in subsonic 
flight . 

If the left side of equation (19) is integrated by parts, the 
general relation becomes 

• -//I“' dxat '//S' ,d " dt (20) 

If the motion of a three-dimensional wing is to be studied, equations 
( 19 ) and (20) must be modified to include an integration with respect 
to y. 

Two further, changes in equation (20) serve to simplify applications. 
In the first place, a symmetry is restored to the expression if two 
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distinct systems of axes are used as in equation (l4); in the second 
place, the pressure relation 


P-Po = 


dcp 

-Po — = 

at« 


aq> 

-Poao r - 

dt 


( 21 ) 


where p Q denotes undisturbed pressure, permits the introduction of 
pressure p in the integrands. The final expression for the three- 
dimensional case is, therefore. 


fff Pi(x 1 ^yi»t 1 )w 2 ( x 1 ,y 1 ,t 1 )dx 1 dy 1 dt 1 = 

Iff P a ( x 2 »y 2 7 't 2 )w 1 (x 2 ,y 2 ,t 2 )dx 2 dy 2 dt 2 


( 22 ) 


where the two motions now follow the same path in reverse directions but 
are referred to the two sets of oppositely oriented axes satisfying the 
relations 


xi = -x 2 + I, 


yi = -y2 + n. 


Zi = z 2 , 


tl = -ts + T (23) 


where t, tj, and t fix the relative positions of the two origins. 

accompanying sketch indicates one 

possible orientation of the axes. z, 


The 


Equation (22) reduces to a much 
simpler form provided further 
restrictions are imposed on the 
upwash functions w 2 (xi,yi,ti) 
and wi(x 2 ,y 2 ,t 2 ) . In order to fix 
the idea, consider the case in which 
the wings have traveled a time 
T* , (T' = T/a 0 ) and a distance 
U 0 T' = MqT. Let the two systems of 

coordinate axes be placed such that % = 0 sets the starting time of 
the forward motion and t 2 = 0 sets the starting time of the reverse 
motion; the two origins are furthermore oriented such that they are at 
opposite ends of the root chord of the common plan form. Equation (22) 
then becomes 



f dt iff Pi( x i»yi^ t i) w 2( x i^yi> t i) dx i d yi = 

0 P l(t!) 

/ m _ _ 

dt 2 f I P 2 ( x 2 ,y 2 , t 2 )wi( X 2 ,y 2 ,t 2 )d x 2 dy 2 
U P 2 (t 2 ) 
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where the functions w x and w 2 have an implicit dependence upon T.. 
If w x and w 2 remain constant for x x + M 0 t 1 = const, or 
x 2 + M 0 t 2 = const., the egression 



x 2 ,y 2 ,T)w 1 (x 2 ,y 2 ,T)dx 2 dy 2 (2i i) 


follows after taking a derivative with respect to T of the original 
equality. Equations (24) and (l4) are now equivalent in form, with T 
taking the role of an auxiliary parameter. In this way, certain classes 
of- unsteady motions can he treated simultaneously with steady motions. 

In the applications to follow it will be convenient to introduce 
into equation (24) upwash functions of the indicial type; that is, 
functions that are zero up to a fixed time and, after experiencing a 
finite discontinuity, remain constant for all subsequent values of time. 
Such indicial or step variations can be assumed, say, for angle of attack 
rate of pitch, and rate of roll since they satisfy the requirements 
underlying the derivation of equation (24). This choice of functions 
will prove to be advantageous in that the integrals of the responsive 
pressures will yield results relating the wing characteristics. Theorems 
to be given later will speak specifically of steady and indicial motions. 
It is to be understood, however, - that the indicial results can be further 
extended when the same wing is assumed to be executing identical motions 
in forward and reverse flight. Thus, by means of Duhamel's integral 
(see, e.g., reference 19 ), if f(t) is the response in the wing character 
istic to a step variation in w at time t = 0 , the response to an 
arbitrary variation with time of w can be written 

F(t) = ^ J' f(t - T)w(T)dr (25) 

0 

If it is known, for example, that the lift per unit angle of attack is 
the same at corresponding values of time for a wing experiencing an 
indicial angle-of -attack change in forward and reverse flight, it fol- 
lows that the build-up of lift is the same at corresponding values of 
time for all forward and reverse motions, provided the time histories of 
the motions are the same. The equivalence of lift would thus be estab- 
lished, for instance, for oscillatory variations ■ in angle of. attack. 


APPLICATIONS 


The results of the foregoing analysis may be employed to determine 
a number of special theorems that are particularly useful in the calcu- 
. lation of the aerodynamic characteristics of twisted wings and of wings 
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in nonuniform downwash fields. The theorems apply equally to wings act- 
ing either alone or. ,in certain cases, in combination with other wings or 
with cylindrical bodies having their generators alined with the x axis. 
Moreover, they apply not only to wings in steady motion but also to wings 
performing unsteady motions of the indicial type, or unsteady motions 
derivable therefrom. For wings in more complex unsteady motions, however, 
it will be necessary to refer to the more general results of equation (22) . 
Some problems of this nature will be described at the end of this section. 

The applications to be included are exact within the framework of 
linear theory and involve no further restrictions on the wing plan form 
or Mach number except in certain indicated cases where it will be con- 
venient to use results based on slender wing theory. The examples are 
intended to be representative in nature. 


REVERSAL THEOREMS - STEADY AND INDICIAL MOTIONS 


Reversal theorems are defined here as relations between the aero- 
dynamic characteristics of identical wings executing the same type of 
motions in forward and reverse flight. The results presented in this 
section apply not only to single wings in steady motion but also to com- 
binations of wings, as in cascades or multiplanes, performing either 
steady motion or motions of the indicial type. 


Drag of Symmetrical Nonlifting Wings 


The drag of a symmetrical, sharp-edged wing in linear theory may be 
determined by integrating over the plan form, the product of the pressure 
and the slope in the x direction of the wing surface; when the wing has 
blunt edges with slopes having square-root singularities, these singular- 
ities yield an added contribution (reference l8) . In general, therefore, 
the drag D of a symmetrical section is given by 


D = D e + 2 


//’ (it) 

P- N 'XL 


dS 


( 26 ) 


where D e is the drag attributable to the edges. 

If the subscripts 1 and 2 refer to the same wing in forward and 
reverse flow, respectively, and with the two systems of axes introduced 
in equation ( 23 ), local slopes are related. as follows 


dz i( x i>yi> T ) 


dx. 


dz 2 (x 2 ,y 2 , T ) 


u 


dx. 


( 27 ) 


u 
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Equations (26) and (27), together with the reciprocal relation ( 2k ), 
yield 



Since the geometry of the wing is fixed, the edge contributions are the 
same, 

(D e ) 1 « (D e ) 2 (29) 

and, consequently , 

Dx = D 2 (30) 

which confirms the relation stated in reference (9) • 

\ 

THEOREM: The pressure drag in steady or indie ial motion of sym- 
metrical nonlifting wings is the same in forward and reverse flight. 

Lift on Flat -Plate Wings 


The lift L of a wing may he determined by integrating the differ- 
ential pressure Ap = p i - p u over the wing plan form, thus 


-L 


Ap dS 


(3D 


For flat-plate wings, the local angle of attack of the wing surface' is a 
constant 


a 1 (x 1 ,y 1 ,T) = = const., a 2 (x 2 ,y 2 ,T) = a 2 = const. (32) 

Application of equations (31), (32), and {2k) yields the following: 

'^Ap 1 a 2 dS 1 = f F APgO^dSg = 


La = 
1 2 


or 


r i. 


Li/ a l = L 2 / a 2 


(33) 
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THEOREM: The lift per unit angle of attack of flat -plate wings 
in steady or indicial motion is the same in forward and reverse 
flight . 

This theorem generalizes the relation previously given by Brown (refer- 
ence 5) for steady motion. 


Damping in Roll of Flat -Plate Wings 


The rolling moment L' exerted on a wing, following the usual sign 
convention, is given by 

L ' = -f /yApdS (34) 


The local angle of attack due to rotation about the x axis is 


a i = 


Px’71 


u„ 


a 2 = 


V 2 , y 1 




u 


(35) 


where p' is the angular velocity of roll, assumed constant. Applica- 
tion of equations (34), (35 ) > and (24) yields the following: 

”\ff y ^ ldSl ' E fci7 ysfiPldSl - 




yiAp 2 dS 2 = - 


V 

Uo 



y 2 Ap 2 dS 2 = 


V V 


or 


V/Pi' - V/Pa' 


(36) 


THEOREM: The rolling moment per unit angular rolling velocity 
of flat-plate wings in steady or indicial motion is the same 
in forward and reverse flight. 


Damping in Pitch of Flat -Plate Wings 


Consider a wing, first, in forward flight and pitching with a uni- 
form angular velocity q x about a lateral axis; second, in reverse 
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flight and pitching with angular velocity q 2 about another lateral 
axis. Place each wing in a coordinate system such that the y axis 
coincides with the axis of rotation and designate the distances to the 

moment axes by Xq with proper 

z, q * subscripts as shown in the adjacent 

u f sketch. In such a coordinate sys- 

>» ! tern, the pitching moment Mxq 

exerted on a wing, following the 
Jl!iPA> usual sign convention, is 


(x-x Q ) ApdS 


The local angle -of -attack distributions due to rotations about the 
y axis are 


o-i = — xi, 

u o 


^ ( t \ 
0-2 = 7T x 2 = — U-xJ 

u o u 0 


Application of- equations ( 38 ) > (39)# and (24) yields the following: 


~ ( M i)x 01 + (6 -Xqi)Li = ^ff ( x oi“ x i)'+ (^" x oi) APidSi = 


(l-x 1 )Ap 1 dS 1 = -2 / / x 2 APidSj. = 


— / / XiAp 2 dS 2 


( M 2)xo 3 + ^- X 02) L « 


(Ml) X 0 l + ^O^ L 1 (M 2 } X op + ^- X 0 2 ) K 


This equation indicates that the pitching moment due to pitching veloc- 
ity is, in general, not the same for wings in forward and reverse flight. 
However, if X 01 = X 02 = |, the pitching moment per unit angular pitching 
velocity of flat-plate wings in steady or indicial motion is invariant. 
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SPECIAL RECIPROCAL THEOREMS AND APPLICATIONS 


In the following section, several special reciprocal theorems will 
be derived and applications will be illustrated. Reciprocal theorems, 
in contrast to reversal theorems treated in the preceding section, are 
defined here as relations between the aerodynamic properties of wings in 
forward and reverse flight that have dissimilar camber, twist, and thick- 
ness distributions but have the same plan forms. The motions may or may 
not be similar, although it is assumed in this section that both wings 
are in either steady motion or unsteady motion of the indicial type. As 
noted in the preceding section, the results apply equally to wings act- 
ing alone or in combination. 


Symmetric Norilifting Wings - Steady Motion 


- The problems of paramount interest in the application of the general 
relations are found from considerations of pressure integrals over lift- 
ing surfaces; such problems will be given detailed treatment later. In 
the present section, a brief indication is given of the manner in which 
useful results -can be derived for symmetric wings. The discussion will 
be limited to steady-state, two-dimensional, subsonic pressure fields 
although fairly obvious extensions can be carried out. 

If the geometry of a real symmetric airfoil is prescribed, the 
theoretical pressure distribution exists and is unique. If, however, 
the pressure distribution is prescribed, a real airfoil does not neces- 
sarily exist, but by means of reciprocal relations 'it is possible to 
derive certain conditions of compatibility that need to be imposed. Con- 
sider, therefore, the two subsonic solutions 

u i( x i) =0, Wj/xJ = Uo/ya 2 ^ 2 ; -a<x x <a 

and 

•u 2 (x 2 ), w 2 (x a ) . 

The first solution has square -root singularities in w at each end of 
the airfoil and, correspondingly, singularities in u occur just ahead 
of the point x = -a and just behind x = a. On the other hand, equa- 
tion ( 14 ) certainly applies if w 2 is zero at x = ±a. If the origins 
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of the two systems of sixes are at the same position, it follows from 
equation (l4) that u 2 (x) must satisfy the relation 

(ko) 

This result is useful in the calculation of airfoil shapes involving a 
change in pressure distribution from that of a known reference profile. 
The restriction on w 2 at the nose and tail implies that the derived 
and reference profiles have the same slope and radius of curvature at 
those points. The restriction on u 2 , as given in equation (4o), can be 
interpreted as a condition that must exist by virtue of the fact that 
the drag of an airfoil in two-dimensional potential flow is zero. 

As a second example, consider the solutions 


r a u 

0 -4 - 


,(x 2 )dx £ 


*/a 2 -} 


( \ — tu ° 
u l( x l) ~ “ y 
2a 


v l( x i) 


-txiU Q 
2a «/a 2 -x 1 2 


-a<x 1 < a 


that represent velocity and slope of a thin ellipse of thickness t and 
chord 2a. If W 2 is chosen as above, such that it vanishes at the nose 
and tail of the airfoil, if u 2 is the corresponding velocity distribu- 
tion, and if the two sets of axes are as before, equation (l4) yields 



From this result, together with the general closure condition 






0 


(4l) 


a necessary condition for the closure of the second airfoil is 



Ui( Xg) X£dX£ 

- 0 

J a 2 *-x 2 2 


(k2) 
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As a final example, consider the solutions for -a < xi < a 

Ui( Xl ) = (a-2 Xl ), 

3 v3a 2 


wi(xi) = 


2tU 0 2x i 2 ~ ~ a 2 

3 a/ 3 a 2 ,/a 2 - x x 2 


representing velocity and slope of a thin Joukowsky type airfoil. In 
this case, w x vanishes at the tail and the downwash distribution w 2 
for the reverse wing may have a square -root singularity at the nose. 

The nose of the first wing is, however, blunt and for equation (l4) to 
apply the second wing must have a cusped tail. Under these Conditions, 
equation (l4) yields 


2t up 2 . r a ( 
3 JJ a 2 J~ a 


a-2x a ) 




2xi^ - axx - a 2 

*/a 2 - Xx 2 


u 2 (x 2 )dx 2 


Making the substitution Xx ='-x 2 in this equation and integrating the 
left side by parts, one has 1 


U„ 



2x a 2 + ax 2 
J a 2 -X2 2 


u 2 (x 2 )dx 2 


For all real airfoils with cusped trailing edges, therefore, the area As 
can be expressed as 


As 



2x 2 2 + ax 2 - a 2 
J a 2 - x 2 2 


Us(x 2 ) 

Uo 


(^3) 


Lift - Steady and Indicial Motion 


The reciprocal theorems offer considerable advantage in the calcu- 
lation of the lift of wings having a nonuniform angle -of -attack distri- 
bution or of wings in a stream having nonuniform flow directions. For 
these applications, it is convenient to consider a special form of the 
reciprocal theorem which relates the lift on a wing having arbitrary 
distribution of local angle of attack to that of the flat -plate wing of 
identical plan form in flight in the reverse direction. Since the 
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solution of this latter problem is often known or can be found rela- 
tively easily, the solution of the original problem is facilitated in 
many instances. 

Lift of arbitrarily cambered wings .- Consider two wings of identical 

plan form in flight in opposite 
directions as shown in the adja- 
cent sketch. Wing 1 is arbitrarily 
cambered and twisted and wing 2 is 
flat . 

a i = °a (xi,yi,T), a 2 = const. 

m 

Application of equations (lj-4) and ( 2k ) yields the following: 




THEOREM: The lift in steady or indicial motion of a wing having 
arbitrary twist and camber is equal to the integral over the 
plan form of the product of the local "angle of attack and the 
loading per unit angle of attack at the corresponding point of 
a flat-plate wing of identical plan form in flight in the 
reverse direction. 


Equation ( 45) may be used to derive Munk's integral formula for the 



lift of an arbitrarily curved air- 
foil in subsonic flow. Consider 
airfoils 1 and 2 placed in their 
respective coordinate systems, as 
indicated in the adjacent sketch. 

The angle -of -attack distributions on 
the two airfoils are given by 


dz x 

ai , <x 0 = const. 

dx x 1 ^ 


( 1 * 6 ) 


The loading per unit angle of attack on airfoil 2 is 


Ap 2 _ *KLo /a-x 2 _ 4q 0 / a~+~x7 
a 2 P V a+x 2 3 v a - x 1 


(V7) 
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"1 2 

where is free-stream dynamic pressure ( -A p 0 U 0 ). Substitution 

into equation (47) yields the lift formula ^ 


Ll , - ISO r* /a^ 

P J dxiV a - xx' 

— £L “ 


(48) 


The corresponding formula for the lift of a tier of curved airfoils 
may also be derived similarly from the expression for the loading on an 
equivalent tier of flat airfoils. 

The load distribution per unit angle of attack for a two-dimensional 
supersonic wing is 


Ap a _ 4qo 

a 2 p 


and, from equation (47) > the lift is 


Li = - 


4q Q r a dzx 

g J _ a dx x 


dx 1 


(1*9) 


(50) 


The extension of this result to include supersonic -edged wings with 
straight trailing edges leads to a result given originally by Lagerstrom 


and Van Dyke (reference 20). If, as 
in the sketch, the sweep angle of the 
straight edge is A, the load distri- 
bution per unit angle of attack of 
the reversed wing is 

Ap g _ 4q Q cos A (51) 

0,2 ^I-Mq 2 cos 2 A 



and the lift is 

L i = = f f ~ dxxdyx 

sec 2 A-Mq 2 J J dxx 
Pi 


(52) 


A less obvious application yields the build-up of lift with time 
of an arbitrarily cambered two-dimensional supersonic airfoil starting 
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impulsively from rest at a constant speed. Rewriting equation ( 45) , 
lift is 


* • -L 


aa-MoT dz-^x^T) Ap 2 (x 2 ,T) 


-MqT 


dxi 


a 2 


dx 2 


(53) 



The sketches show the final 
positions of the airfoils 
relative to each other for various 
values of T. In reference (21) 
the expressions for the loading in 
equation (53) are given. Over 
the intervals denoted by 
©, ©, and © in the sketch, 
these expressions are 


Region ©) 
Region ( 2 ) 


Ap 2 _ 4q 0 • 
a 2 Mq 

H- 


Ap 




«a/ m q 2 -1 


arc cos 


M q x+T 

x+MqT 


S 


Mp - 1 
Mo 


(i 


+ arc sin - 
T 


Region © 


4Pa = 

a 2 J Mq 2 -! 


J 


Lift on a wing in a nonuniform downwash field .- The reciprocal 
theorem of equation ( 45) can also provide a particularly good method of 
determining the lift on a wing in certain nonuniform downwash fields of 
known structure. Such problems arise whenever a wing acts in the 
presence of other wings, bodies, or propellers but is always of prime 
concern in the determination of the lift on a tail acting in the down- 
wash field of a wing. In most problems, the downwash velocities at the 
position of the tail may be considered to be constant in the longitudi- 
nal direction and to vary in the spanwise direction, thus 


aa = My^T), 


02 = const 


(55) 
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and 



where Z 2 is the span load distribution associated with the load distri- 
bution Ap 2 . Summarizing, the lift in steady or indicial motion of a 
wing in a downwash field which varies across the span is equal to the 
integral over the span of the product of the local angle of attack and 
the span loading per unit angle of attack at the corresponding spanwise 
station of a flat-plate wing of identical plan form in flight in the 
reverse direction. This statement generalizes the result given recently 
by Alden and Schindel (reference 22) for steady flow about wings having 
supersonic leading and trailing edges and streamwise side edges. 

As an example, consider the problem of determining the lift on a 
wing at a geometrical angle of attack of zero resulting from the pres- 
ence of an infinite line vortex of 
strength r extending in the flight 
direction. The wing will be con- 
sidered to have such a plan form that 
its span loading when in flight in 
the reverse direction is elliptic. 

The notation is as shown in the 
adjacent sketch. For this problem, 
therefore, the span loading of the 
wing in reverse flight is given by 




The local angle of attack of the original wing due to the presence of the 
vortex is given by. 


r yi-P 

2k u Q f(y r p) 2 + 7 2 ] 


( 58 ) 


Substitution of equations ( 57 ) and ( 58 ) into equation ( 56 ) yields the 
following formula for the lift: 



The lift oh a wing in the vicinity of a number of such vortices may be 
found by superposition. 
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Lift due to deflection of a portion of the wing surface ,- Let a 
portion P* of the surface of wing 1 be deflected a constant angle 6 
and the remainder of the wing be a flat plate alined with the free -stream 
direction. Let wing 2 be 'a flat-plate wing inclined at an angle of 
attack a, thus ( 

5 on P 1 

a,! = < a 2 = const. (60) 

0 elsewhere 

Substitution of equation (60) into (45) yields the following result: 



The lift in steady or indicial motion per unit angular deflection 
of a portion of the wing surface is thus equal to the lift per unit 
angle of attack on the corresponding portion of a flat-plate wing in 
flight in the reverse direction. This generalizes a result given 
previously by Morikawa and Puckett (reference 23) for steady flow about 
low-aspect-ratio wings. 

This rule is very useful in the determination of the lift resulting 
from the deflection of a flap or control surface. This is particularly 
true for supersonic speeds since the loading on the related flat wing is 
often a constant over a large portion of. the area of integration. 

As a further example, consider the case of a low-aspect-ratio wing 
having a straight, trailing edge and mounted on an infinite cylindrical 
body of revolution. The entire wing -body combination is at zero angle 
of attack except for the flaps on the rear of the wing that are deflected 
an angle &. The problem is to determine the lift on the entire wing- 
body combination due to the deflection of the flaps. Slender-wing- 

theory results of reference 2k are 
to be used. The notation is indi- 
cated in the adjacent sketch. 

The solution of this problem 
is particularly facilitated by the 
fact that slender-wing theory indi- 
cates that the loading on wing 2 is 
concentrated on the leading edge as 
shown in the sketch. Therefore, the 
lift of wing 1 is found by integrating the span loading curve of wing 2 
over the portion of the span between y = f and y = s. Thus 
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The span loading on wing 2 is given by 



12 - 

a 2 


a2 




i 2 (y a ) = “ p 

The lift due to the deflection of the flap is therefore 


L flap “ 4( lo 5s2 

I C 1 + ?) 


C 1 -fij c - /i 1 - p) (i 


a 4 , 

B* > + 


. 1 - 2(f/s) 2 + (a/s) 4 

arc sin , , , / + 

1 - (a/s) 4 


a 2 

— arc sin 


[l + (a/s) 4 ] (f/s) 2 - 2(a/s) 4 
(f/s) 2 [l-(a/s) 4 ] 


(63) 


(64) 


A plot of the results is shown in the 
adjacent sketch. The lift per unit 
angle of flap deflection (l/ 8.)^ 
has been nondimens ionali zed by 
dividing by the 'lift per unit angle 
of attack (L/a)^_g of a slender 
wing-body combination of identical 
plan form,. From reference 2k, 

(l/o 0 w _ b is given by 


(L/a,)y_B - 2itq. 0 s 


(65) 



Rolling Moment - Steady and Indicial Motion 


The calculation of the rolling -moment characteristics of wings 
having a nonuniform angle -of -attack distribution or of wings in a stream 
with nonuniform-flow directions can be performed in many cases through 
use of the reciprocal theorem in a manner analogous to that described 
for the lift characteristics in the preceding . section. In every case, 
the rolling moment of the given wing will be related to the lift on a 
rolling flat -plate wing. For the sake of simplicity, all of the present 
examples will be confined to the case where the rolling moments are 
^v&luated about the x axis, considered to lie in the plane of symmetry. 
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Rolling moment of arbitrarily cambered wings.- Consider two wings 


** 



of identical plan form in flight in 
opposite directions as shown in the 
adjacent sketch. If the local 
angle -of -attack distribution of 
wing 1 is arbitrary and that of 
wing 2 varies linearly with . y 
(which might be likened to either 
a wing with linear twist distribu- 
tion or to a flat-plate wing 
rolling about the x axis with 
constant angular velocity P 2 *)> 


ai = a x U lt y x ,T), 



and the following relations can be written: 


( 66 ) 



THEOREM: The rolling moment in steady or indicial motion of a 
wing having arbitrary twist and camber is equal to the integral 
over the plan form of the product of the local angle of attack 
and the loading per unit (p 2 , /u 0 ) at the corresponding point- of 
a rolling flat -plate wing of identical plan form in flight in 
the reverse direction. ✓ 

Applications of this theorem follow in a manner very similar to that 
described previously for the corresponding theorem regarding lift. 

Rolling moment on a wing in a nonuniform downwash field .- Consider 
a wing placed in a flow field in which the downwash velocities at the 
position of the wing are constant in the longitudinal direction and vary 
in the spanwise direction. The related wing is again a flat-plate wing 
rolling with an angular velocity p 2 f . 
section, thus 


as described in the preceding 


oc 2 — 


P ^2 


a l = a-x (Yi,T)> 


( 68 ) 
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The rolling moment of the first wing is then given by 



(69) 


or > -*- n words, the rolling moment in steady or indicial motion of a wing 
in a downwash field which varies across the span is equal to the integral 
over the span of the product of the local angle of attack and the span 
loading per unit (p 2 '/ u o) at 1 the corresponding spanwise station of a 
^^Lling flat -plate wing of identical plan form in flight in the reverse 
direction. • 


^^Lling moment due to deflection of a portion of the wing surface » — 
Let a portion P* of the surface of a wing be deflected a constant * 
angle 5 and the remainder of the wing be a flat plate alined with the 
free-stream direction. The related wing is a flat-plate wing rolling 
with angular velocity p 2 ' 


oa 


[ 8 on P* 
j 0 elsewhere 



(70) 


Substitution from equation (70) into ( 67 ) yields the following result: 



(71) 


Thus, the rolling moment in steady or indicial motion due to a 
given angular deflection of a portion of the wing surface is equal to 
the lift per unit (p 2 '/lJo) on the corresponding portion of a rolling 
flat-plate wing of identical plan form in flight in the reverse direc- 
tion. 


As an example, consider a wing-body combination consisting of a 
low-aspect-ratio wing having a straight trailing edge mounted on an 


infinite circular cylinder as shown 
in the adjacent sketch. The body is 
at zero angle of attack, the right 
wing P' is deflected an angle 8 
and the left wing P' ' is deflected 
-8. The problem is to determine by 
means of slender-wing theory the 
rolling moment exerted on the entire 
wing -body combination. The notation 
is indicated in the adjacent sketch. 


V 
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a-i = 


- 5 on left wing, P* ' 
0 on body 

+ 5 on right wing, P* 


a 2 = 


Pa'yg 


(72) 


Since slender-wing theory indicates that the loading on wing 2 is con- 
centrated on the leading edge, the rolling moment of wing 1 can he found 
by integrating the span loading on wing 2. 


p 2 ' 


AP g \ 
2 ’/Uo ) 


,-a 


dS 2 = - & 


A) dya + 6 ^(^fe) 


dyo = 


28 ' C-i" ] ^ 


where, from reference 25 , 

Zo 


Pp’/Uo 


= lo s 


n 2 2as 

1 + 7 arc cos 55^5 


) ( y + r) / s2 + S - ** - % 


2 / a 2 \ 2 . (y 2 + a 2 )(s 2 - a 2 ) 

= ( y - -y ) arc cosh 

\ J (y 2 - a 2 )(s 2 + a2) 


(73) 


The resulting expression for the rolling moment is 


8 E 


E(cp,k) 






4 ' ' 2 _ _ 2P \ . 2 _r>2 


— arc cos 


arc cos — -*■ i - — 

1+Fp J 3n 

— R(l-R 2 ) 


”)-fH 


(l-R 2 ) + 


l+R* 


0 + 5 r(i - 


2R 


— In - -g + — (1 - R 2 ) (l + \ arc cos 
3* 1 + R 2 3 V * 1 + R 


where 
R = a/s 


( 74 ) 


9 = arc sin 


Vl + R 2 


k = «/l - R 4 
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A plot of the results is shown 
in the adjacent sketch. The rolling 
moment has been nondimensionalized 
by dividing by the value correspond- 
ing to that of the wing alone 
(R = 0). 


Pitching Moment - Steady 
and Indicial Motion 


A number of useful relations 
regarding the pitching -moment charac- 
teristics of wings may be found _ 
through application of the reciprocal 
theorem. Since the general procedure 
is closely analogous to that of the 
preceding sections, the following 
discussion will be brief. 
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Pitching moment of arbitrarily . cambered wing«- Consider the problem 


of determining the pitching moment 
Mi about the origin of wing 1 pos- 
sessing an arbitrary distribution of 
camber. The related wing in flight 
in the reverse direction, -wing 2 , is 
a flat -plate wing of identical plan 
form pitching about the moment axis 
of wing 1 as indicated in the sketch, 
thus 



= a x (xx,yx,T), 


a 2 


q.^1 q2(x 2 -|) 

“Uo ~o 


The pitching moment of wing 1 is given by 



( 75 ) 

( 76 ) 


THEOREM: The pitching moment in steady or indicial motion of a 
wing having arbitrary twist and camber is equal to the integral 
over the plan form of the product of the local angle of attack 
and the loading per unit (q 2 /U 0 ) at the corresponding point of 
a flat -plate wing of identical plan form in flight in the 
reverse direction and pitching about the moment axis of the 
first wing. 
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The necessity for pitching wing 2 about the moment axis of wing 1 
may be removed by considering wing 2 to be re-expressed in terms of two 
component wings having angle -of -attack distributions given by 

0-2 * ( X 02 1 ) 

a s t ' = — = const. (77) 

U Q 

Wing 2' is thus pitching with angular velocity about an axis at 

x 2 t = x 02 i and wing 2’ ' is a flat -plate wing at a constant angle of 
attack. The pitching moment on wing 1 is then given by 


a. 


_ V 


( X 2. - 


v 02' 


U r 



Applications of pitching -moment theorem. - Since the application of 
equation (76) or (78) to problems analogous to those discussed in the 
preceding sections proceeds in a straightforward manner, the following 
discussion will be abbreviated. 


If wing 1 is a flat-plate wing, is constant and equation (78) 
simplifies to 


Ml 

a ! 


<Ls' /Uo + 


(x o2 » - I) 


Lgi i 

a 2 r t 


(79) 


where L 2 ' is the lift on wing 2' pitching about x 2 1 = Xq 2 1, and L 2 m 
is the lift on an inclined flat-plate wing. Equation (79) may be 
expressed in terms of conventional stability derivatives as follows: 


( Cm a\ “ ( C L a ) 2 * + 



I 1 


(80) 
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An application of this result to 
unsteady flow problems is indicated 
in the adjacent sketch obtained from 
indie ial- lift and pitching-moment 
results of reference 21 . This 
sketch shows the growth of lift and 
pitching moment on triangular wings 
with supersonic edges at a Mach 
number of 2 following indie ial angle - 
of -attack and pitching-velocity 
changes. In these results/ the 
rotation and moment axes are always 
at the leading edge or apex, there- 
fore Xqs 1 = 0 and £ = c Q . It may 
be seen that the three curves are 
related in the simple linear manner 
indicated by equation (80) . 

If is independent of x and 

varies only in the spanwise direction, 
that is, if <x x = a 1 (y 1 ), the pitch- 
ing moment on wing 1 is given by the 



following equation, analogous to equation (56) for lift. 



If a portion P f of the surface of wing 1 is deflected a constant 
angle 6. and the remainder of the wing is a flat plate alined with the 
free -stream direction, the following relations hold: 

6 on P T . 

0 elsewhere 

and 
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RECIPROCAL RELATIONS INVOLVING MOTION INTO A GUST 


All previous applications that have been considered were derived 
from equation (24). In the present section, the more general equation 
(22) will be used to develop two theorems which relate the build-up of 
lift on a wing entering a gust and the build-up of circulation on the 
same wing moving indiciall'y but in the opposite direction. The rela- 
tions to be obtained hold for the Mach number range for which the wave 
equation applies. Under the special assumptions of incompressible flow, 
the results in two dimensions establish a direct connection between the 
circulation function calculated by Wagner (reference 26) and the gust 
lift-curve calculated by Kussner (reference 27). A proof of the agree- 
ment between these latter results has been given by Sears (reference 28). 


Two-Dimensional Flow 


A flat plate is assumed to be moving in two modes of motion: In 

the motion associated with the axes (x a ,z 2 ,t 2 ) the wing starts at time 
zero (t 2 = 0) and moves at a constant velocity U Q and at a constant 
angle of attack; the motion associated with the axes (x 1 ,z 1 ,t 1 ) starts 
at time zero (t x = 0) with the wing moving in the opposite direction at 
a velocity U 0 and entering a sharp-edged gust. The gust exists for 
all xi less than zero and has a vertical velocity Wg = - a,g U 0 . The 
two wings, therefore, have angles of attack as follows: 

a 2 (x 2 ,t 2 ) = a 2 = const. for -M 0 t 2 < x 2 < 2a-M 0 t 2 , t 2 >0 


a^Xi,^) = a,g = const. 


for 


0 >x z >-M 0 t x , . 0 <t x ^ 2a/M 0 

2a-M 0 t! >Xx >-Moti, 2a/M 0 ^ t x 


The two-dimensional form of equation (22) yields 



■¥1 


2a " M o t i 


Ap g (x 1 ,t 1 )dx 1 = P 0 a o a g J f A <*S 

A 2 


(83) 


where the region A is bounded by the lines xi = 0, x x = 2a - Mq^, 
ti = T, and x-l = -M 0 t 1 . The integral on the right can be rewritten as 
a line integral by means of the identity 



cos(t,n) 


ds 
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and equation (83) becomes 


l 2 / L g (ti)dti = p 0 a 0 a g / A 9 2 cos (t 2 ,n)ds 


•where the line integral extends around the boundary of the region A. 
Since A92 vanishes on the lines x 2 = -M 0 t 2 and t 2 = 0 , the equation 
becomes 


a 2 / 'L g (ti)dt 1 = p 0 a 0 a„ 


A 9 2 ( x 2 . 


2a-M 0 T 


\ ' 


Differentiation with respect to T yields 


a 2 L g (T) = P 0 U 0 agACp 2 ( 2 a-MoT,T) 


The discontinuity in 92 is evaluated at the trailing edge at time T 
and is therefore equal to the circulation ^2 of the airfoil at time T. 
The equality thus becomes 

L g (T) t 2 (t) 

T = p o U o — “ ( 84 ) 


ob„ P ° U ° a. 


THEOREM: The circulation per unit angle of attack of a flat 
plate moving indicially with a velocity U Q is proportional 
to the lift per unit <Xg of the plate entering a sharp-edged • 
gust having a uniform vertical velocity equal to Wg = -a,g U Q * 



0 5 10 15 0 5 10 15 0 5 10 15 

Half -chords traveled Half-chords traveled Half-chords traveled 

M o s 0 M c - 0.8 M 0 s 146 


In the accompanying sketches, the time variation of these variables, 
as well as the lift of the indicial wing, is indicated for low speed and 
for flight Mach numbers equal to 0.8 and 1.46 as determined from refer- 
ences 21, 26, 27, and 29* 
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Three-Dimensional Flow 


The extension of the above results to three dimensions is 
straightforward. The origin of the xj., yi, Zj., ti axes is assumed 
to be initially at the foremost point of the wing in reverse motion. 
The two wings have, respectively, angles of attack a 2 = const, over 
the reverse moving plan form for all values of time and a.g = const, 
over the region occupied simultaneously by the forward moving wing and 
the gust. Equation (22) gives 


I dt x If APg U^y^ti) dx x d yi = 

0 p(tl) 

r r r dv 2 (x 2 ,y s, t 2 ) 

Po a o a g JJJ A- 


dt . 


dx 2 dy 2 dt 2 


(85) 


The integral on the right can be rewritten as a two-dimensional surface 
integral by means of the identity 

,'p cos (t,n)dS dxdydt 

and equation ( 85 ) becomes 

a 2 J L g (*1) dt i = Wg/ ^ A c P 2 cos ( t 2 ^ n ) dS 2 


where the integral on the right extends over the boundary of the volume 
in x 2 , y 2 , t 2 space occupied by the wing and the gust. The value of 
AT 2 must, of course, vanish on the leading edge of the wing and at 
ta = 0. In order to fix the limits of integration, suppose 'the wing is 
symmetrical about its longitudinal axis and let the leading edge of the 
forward wing be given by the equation 

x i = f (±y x ) - M 0 t 1 or y x = ±s(x x + Mo^) 

For the reverse wing and its coordinate system, this edge, which is now 
the trailing edge, is 


X 2 = c r - Mot 2 - f(±y 2 ) 


or 


,y 2 = ^(c-r - x 2 - M D t 2 ) 
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where c r is the root chord. The reverse-flow integrals then become 


*2>y Si 


c r “*2 

Mq 


T c r s(c r -x 2 ) 

a sf L g (t 1 )dti = p 0 a 0 ag /. dx 2 A9 2 

O c r -M 0 T -s(c r -x 2 ) 

where s(x 2 ) is the local half -span of the wing. Differentiation with 
respect to T yields 

siMoT) 


<ly 2 


a 2 L g (T) = p 0 U Q a g 


ACpc 


-s(MqT) 


f (±y P )"i 

c r -M 0 T, y 2 , T - — g-S- 


dy 2 (86) 


The discontinuity in cp is thus to be integrated spanwise at the rear- 
most point of the indicial wing; this follows from the relation 


Acp. 


x,y, t - 


f(±y 2 ) 

Mo J 


= a <Po (^y^t) j 


T - £ T^< t < T 


which fixe 8 the vorticity in the wake of the wing once it is shed from 
the trailing edge. 


It remains to mention the nature 
of the limits ±s(MoT) . As shown in 
the sketch, the span width of the 
vortex wake at the trailing, edge is, 
during the early stages of the motion, 
dependent on the local span width of 
the wing. The width 2s(MoT) of wake 
is, in fact, equal to the maximum 
width of the portion of the first 
wing that lies within the gust. 

After sufficient time has passed for 
the vortex wake of the indicial wing 
to develop its full span width at the 
trailing edge, s(MoT). becomes s Q or 
semi span of the wing. The nature of 
the wing and wake is indicated in the 
final sketch. 




Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Feb. 19, 1952. 
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